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On 

0^ | There have been many declarations of the death of the Weinberg model of 

spontaneous CP violation. Previous studies, before the recent measurements 
§ H ' of e'/e, indicated that the model could not accommodate the experimental 

values on e in K°-K° mixing, the neutron electric dipole moment (EDM), 

•rH , 

the branching ratio of b — > 57 and the upper limit on e'/e. We point out that 

5-H ' 



these studies were based on optimistic estimates of the uncertainties in the 
calculations and that when more realistic estimates of these errors are used the 
Weinberg model cannot be conclusively ruled out from these considerations 
alone. 

Here we use these realistic error estimates to analyze the present situation 
of the Weinberg model. The latest results from Belle and BaBar on sin 2/3 
allow the small values of this parameter which occur naturally in the Weinberg 
model. However, in this model, the recently measured value of Re(e'/e) = 
(1.92 ± 0.25) x 10~ 3 cannot be made compatible with the branching ratio 
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B(b — > sj) = (3.15 ±0.54) x 10 4 . As a result we conclude that the Weinberg 

model is now confidently and conservatively ruled out. 
ll.30.-j, 12.60.Fr, 12.15.-y 
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I. INTRODUCTION 



The origin of CP violation remains one of the outstanding problems of modern particle 
physics. Although the Standard Model (SM) of CP violation based on the Kobayashi- 
Maskawa (KM) mechanism is consistent |1J] with observations of CP violation in K$ and 
Kl mixing || and in K$,l kit decay amplitudes ||, there are intriguing hints, from 
consideration of baryon asymmetry of the universe ||, that other sources of CP violation 
may exist. These non-standard sources of CP violation could occur as well as, or instead of, 
the SM source. 

Models based on additional Higgs bosons |||6| provide alternatives which explain the ex- 
isting laboratory data and produce the large CP violation required for baryon asymmetry 
Q. Such models also allow CP symmetry to be broken spontaneously |||8| and therefore give 
an interesting explanation of the origin for CP violation. The minimal model of this type 
satisfying the requirement of vanishing tree level flavour changing neutral current ("natu- 
ral flavor conservation" ) is a model of spontaneous CP violation with three Higgs doublets 
proposed by Weinberg |J and refined by Branco H . We shall refer to it as Weinberg model 
from here on. 

It has frequently been claimed that the Weinberg model is in conflict with the data on 
the following: 

• the value of sin 2(3. 

• the ratio e'/e, 

• the e parameter in K° — K° mixing, 

• the neutron electric dipole moment (EDM), 

• the branching ratio of b — > 57, 

First of all, the Weinberg model predicts small values (less than 0.05) for the parameter 
sin 2(3 which here we define as the parameter that characterizes the CP violation in B — ► 
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J/ipKs decays 0,0. The present results |TT[ from Belle (0A5± H(st&t)± Q %l(sys)) and 
BaBar (0.12 ± 0.37 ± 0.09) on sin 2/3 allow such small values of this parameter, although 



the earlier combined result from ALEPH, OPAL and CDF(0.91 ± 0.35) |T2j favoured larger 
values. Therefore these considerations do not rule out the Weinberg model, and to determine 
whether or not the model is consistent with present day data we have to turn to the other 
observables. 

In early discussions of CP violation in the neutral kaon system, it was assumed that only 
the short-distance contributions from the CP violating AS" = 2 box diagrams, due to either 
two charged Higgs particles, or one charged Higgs plus one W exchange, are responsible 
for the measured value of e, then, since the charged Higgs couplings to light fermions are 
proportional to the fermion masses a very large CP violating coupling is required to fit 
the data. If this same CP violating parameter is then used for the calculation of e'/e, the 



contribution is much larger than the experimental value [p^ , ^4| . It was later shown that 
there are important long distance contributions to e due to the Higgs induced As = 1 
CP violating operators, and that if these are taken into account, the model can be made 



consistent with the observed e and e'/e fli"5|,|16 



Previous studies have claimed that although the Weinberg model is consistent with CP 
violation in kaon system, it has problems with the neutron electric dipole moment (EDM) 
T7fl and the branching ratio for b — > S7 |9]Jl0|. These analyses have relied on optimistic 



estimates of the accuracy of the calculations of relevant hadronic matrix elements [|18[ 



After briefly describing the Weinberg model of CP violation, we review the current 
experimental and theoretical status of each of these observables individually, and then discuss 
the constraints placed on the Weinberg model by the total ensemble of data. We find that 
it is possible to rule out this model definitively only if one imposes simultaneously the new 
constraints from the recently measured value of Re(e'/e) = (1.92 ± 0.25) x 10" 3 || and the 
branching ratio B{b -> try) = (3.15 ± 0.54) x 10~ 4 @. 

Therefore, because of the stronger constraints imposed by recent data, one can now 
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declare that the Weinberg model of spontaneous CP violation is ruled out in spite of the 
relatively large hadronic uncertainty. 



II. THE WEINBERG MODEL OF CP VIOLATION 

In the Weinberg model, three Higgs doublets are introduced. The spontaneous break- 
down of gauge symmetry then induces massive Higgs eigenstates for the charged Higgs 
particles, and introduces a mixing matrix specifying the interaction eigenstates of the Higgs 
doublets in term of the mass eigenstates. Being a 3 x 3 mixing matrix between charged parti- 
cles, the matrix contains exactly one irreducible complex phase, thus inducing CP violation. 
After spontaneous symmetry breaking, there are two physical charged and five neutral Higgs 
particles. As CP violation in flavour changing processes in this model is dominated by ex- 
change of charged Higgs particles, we concentrate our attentioin on this contribution. The 
interaction Lagrangian for the coupling of the two charged Higgs (Hf and H}) to fermions 
2Cj can be written as 



L = 2 3 / 4 G 1 J 2 U[V KM M D (a 1 Ht + a 2 H + )R + M v V KM {p x Ht + foH^LjD + H.C. , (1) 

where R{L) = (1 ± 7s)/2, and My^ are the diagonal up and down quark mass matrices. 
The parameters a« and which satisfy lm.(otif3{) = — Im(a2/32), are obtained from the 
diagonalisation of the charged Higgs mass matrix . The KM matrix elements Vij can be 
made all real at tree level as a consequence of spontaneous CP violation. 



III. sin 2(3 IN THE WEINBERG MODEL 

It is well known that CP violation in B decays will eventually provide crucial constraints 
on models of CP violation when sufficient data is available. This is especially true of the 
gold-plated mode B — > J/tpKs- 

In the Weinberg model, CP violating contributions to the decay amplitudes and to 
B° — B° are both proportional to lm(cti/3l). The Higgs contributions to the decay amplitudes 
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are suppressed, relative to the SM contributions, by additional factors of m c mi J /m 2 H , while 
the mixing is supressed by a factor of m 2 /m 2 H . These suppression factors lead to small CP 
violating phases and result in a very small value for sin 2(3 [§,0], | sin 2/?| < 0.05. 



The ALEPH, OPAL and CDF data reported in 1999 gave sin 2(3 = 0.91 ± 0.35 12 



which is in conflict with the above limit at the 2a level. However at ICHEP2000, Belle 



and BaBar reported preliminary results [11], which when averaged with the above, give 



sin 2(3 = 0.49 ± 0.45, consistent with the above limit at the 75% level. Thus the present 
sin 2(3 measurements do not rule out the Weinberg model. 

IV. e'/e IN THE WEINBERG MODEL 

The dominant contribution to e'/e in the Weinberg model is from the flavor changing 
gluonic dipole interaction given by M 



H(sdg) = ig s fm s sa^\ a G^{l - 7e )d, 

f = ^^V ls V ld \m{al^){F^m 2 /m 2 Hi ) - F 3 (mf /m^))^, 

1 X r 3 „ 1 



where i is summed over u, c, t and r\ g = (a s (m^) / 'a s (fi)) 14: ^ 33 ~ 2nf ^ is the one loop QCD 



correction factor |21] in which nj being the number of quark with mass less than /i. To 
obtain this correction factor we will use one loop running for a s with the starting value 
ots{mz) — 0.119. The contribution to e'/e is dominated by the lightest charged Higgs 
exchange. In our later discussions, we will assume is the lighter one and the other is 
very heavy and its effects can be neglected. 

Theoretical analyses for e'/e are conventionally carried out in terms of the isospin am- 
plitude Aj for K — > mi. Expressing e'/e in terms of Aj, one obtains 

where u = ReA 2 /ReA ta 1/22.2. 



The dominant gluon dipole operator Eq (0) generates a non-zero value only for A . 
Calculating the decay amplitudes is our most difficult task, because of our poor under- 
standing of the strong interaction at low energies. Theoretical calculations for the real part 
of the amplitudes can be easily off by a factor of two to three. For this reason we use 
the experimental value for ReAo = 33.3 x 10~ 8 GeV -2 to minimize the error in the cal- 
culation of e'/e. But we still have to calculate Iim4 . This requires the evaluation of the 
matrix element < (ttti)q\0\K >. Here (vr7r) indicates the isospin 1 = component and 
= g a m a 8X a <T ia ,G^(l-j B )d. 

A naive PCAC calculation |2^j23|1 gives 



((inr)o\0\K) = -2^3/2(m 2 m s /(m u + m s )){m 2 K f K / f w ) (4) 



with illy ps lGeV . A bag model calculation [|14[] of 

Ark =< n°\sX a a^G^ (1 - ^)d\K > (5) 

gives 

A Kn ^ 0.4GeV 3 , 



and with the use of current algebra this fl6j gives a value for {{tx-k)q\0\K) similar to that 
of Eq (§. 

It was later realized that the above result is incorrect because an important "tadpole" 
contribution due to the K-vacuum transition caused by the same operator had been ne- 
glected. This contribution cancels the above PCAC result exactly 0,01 • m a chiral pertur- 
bation theory approach, this means that the leading order contribution vanishes as expected 
from the Feinberg-Kabir- Weinberg Theorem ||16| . A non-zero value for < (7T7r)o|C|if > can 



only be generated at p A order in chiral perturbation theory, and can be estimated to be E3 



< (nn) \O\K >= -11 Jl-J^—& mlmlB , (6) 

V 2 m s + m d 

where B is a fudge factor representing the potential uncertainty in the above estimate. We 
assume that B is of order 1. 



Using this matrix element, we obtain 

= 1.7 x 10 7 (GeV 2 )/B , (7) 

where the numerical value follows from the experimental values for e and ReA , and the 
isospin breaking correction factor Vt, q+r] i = 0.25 given in ref p4] 



To produce the recently observed value for e'/e within 3a, fB has to be in the range 
+ (0.69 ~ 1.57) x 10 -10 GeV -2 . For a given Higgs mass, the CP violating parameter Im(ai/3*) 
is determined by the value of /. Only f'B is determined by the data. Remember that the 
different leading order contributions cancel each other. However, numerically the value 
obtained from Eq (|7|) with B = 1 is not much smaller than the individual leading terms 
before cancellations, suggesting that B = 1 is the maximum value of B and hence that the 
corresponding low value for f,f = 0.69 x 10 -10 GeV -2 , represents its probable lower bound. 
Nevertheless, we conservatively allow B to vary from 0.5 to 2 in our estimate to account for 
possible uncertainties p5| . The most conservative range for / is then (0.35 ~ 3.1) x 10 -10 
GeV -2 . Thus / smaller than 0.35 x 10 -10 GeV -2 in magnitude is unlikely to generate e'/e 
as large as observed. Note that / is positive because B , calculated in Ref f23|l , is positive. 

V. e IN THE WEINBERG MODEL 

A successful model for CP violation must to be able to produce the experimental value 
for e. In this model the short distance AS = 2 interaction gives too small a value for e, 
and the dominant contribution actually comes from long distance effects, which in turn are 
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generated by CP violation due to the gluonic dipole interaction. Following Ref. [|Tg 
assume the contribution to e is from ir, rj, rf poles with one CP conserving and one CP 
violating K to 7T, rj, rj' transition. One has [|TBJ 

V2m K Am L _ s (m 2 K - ml) ' 



S 



K = 1 + 



m 



K 



m; 



m 



« — mS! 



rn 



K 



K 



/// 



3(1 + *) cos# + 2W^psin# 



^(1 + 5) sin0-2W|ocos0 



(8) 



where Ami_s is the mass difference of the long and short lived neutral kaons, 9 is the rj — rj' 
mixing angle, and 5 and p parameterize SU(3) and U(3) breaking effects, respectively. In 
the SU(3) limit, 5 = 0; in the U(3) limit, p — 1. H^k is the CP conserving AS = 1, K — n 
transition amplitude which is determined from current algebra to be |L6|| H n x = 2.578 x 10~ 8 
GeV 2 . H-xK — ] '9s m s^-Kn is the CP violating K — ir transition amplitude due to the gluonic 
dipole interaction, expressed in terms of the matrix element A^n of Eq ([5]). Here the QCD 
coupling constant g s is evaluated at the Kaon scale and is not well determined. Following 
Ref. f26fl , we use g s = 47t/a/6 in the matrix element calculation. For m s , one should take 
the values used in conjunction with the models used to calculate relevant matrix elements, 
for example in the bag model calculation it is in the range of 0.3 to 0.5 GeV ||27|| . 

The parameters, 6, 5, p, and the theoretical calculation of Ak-k, introduce uncertainties. 



At present, there are two possible values —11° and —22° [gg] for the mixing angle 6. The 



SU(3) breaking parameter 5 is theoretically estimated p9| . The best fit Kl —> 77 and 
K — > 7T7T7 gives a U(3) breaking parameter p = 0.78. However, p in the range of 0.7 ~ 1.3 
is not ruled out. Using 9 = —22°, 8 = 0.17 and p = 0.78, one obtains k = 0.2. Most of the 
previous calculations used this value for k. However, the value for k is very sensitive to the 
specific values of the parameters involved, for example, with 9 = — 11°(— 22°), p = 1.3, and 
5 = 0.17(0.0), k is approximately -0.9 (-0.95). The magnitude of k can change by a factor 
of four or five. Within the allowed parameter space, k can vary between 0.2 to -1.0. We 
note that the sign of k changes in the allowed range of parameters which implies that the 
relative sign of e and e' can change. 

The uncertainty in the value A wK is also quite large. A bag model calculation gives A wK = 
OAGeV 3 [14]]. Major sources of uncertainty include the determination of the numerical values 
of a s at the kaon decay scale, of the bag radius R and of the strange quark mass m s in bag 
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model |27[]. A factor of two to three times increase in A^k is n °t ruled out. In view of 
these uncertainties, we consider («, A^-k) in the rectangle with corners (0.2, 0.4GeV 3 ) and 
(—1.0, 1.2GeV 3 ) to be allowed by present experimental and theoretical estimates. Of the 
two extreme values are, set a) (0.2,0.4GeV 3 ) and set b) (-1.0, 1.2GeV 3 ), set a) is the mostly 
used one in the literature, while set b) represents the most conservative values for k and 

We find that if the parameters of set a) are used to fit e, the parameter / is determined 
to be 2.56 x 10~ 10 GeV~ 2 . However, if set b) is used, / can be negative and as small as 
0.17 x 10~ 10 GeV -2 in magnitude. There are solutions for e with / in the ranges (0.85 ~ 
2.56) x 10~ 10 GeV" 2 and also with / near -0.17 x 10~ 10 GeV" 2 . The allowed range of / 
associated with e is thus quite large, and it has a large overlap with that determined from 
e'/e for positive /. 

VI. THE NEUTRON ELECTRIC DIPOLE MOMENT 

The experimental bound on the neutron EDM, d n , has been used to provide restrictions 
on the model, and has been claimed to rule it out [|17] . The neutron EDM can be generated 



by the exchange of neutral and charged Higgs particles pIH34]]. It is not impossible that 



these contributions may cancel each other and result in a very small neutron EDM. Here 
we will not entertain this possibility. We will instead single out the variously potentially 
large valence quark contributions and require that each of them satisfies the experimental 
constraints. 

The contribution of charged Higgs exchange to the neutron EDM is well constrained by 
fixing the CP violating parameter Im(ati/3*) to fit e'/e and e. The contributions from neutral 
Higgs exchange are much less constrained. Even in the charged Higgs case we need to use 
the theoretical expression for / to extract Im(oti/3*), and this introduces a sensitivity to the 
values of the KM elements because the internal charm and top contributions are comparable 
and can add constructively or destructively depending on the relative sign of combinations 
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of the KM matrix elements. This also introduces uncertainties in the calculations. The case 
where the contributions tend to cancel will result in a large Im(o! 1( 5j) and lead to difficulties 
with other data as discussed below. We will use values of the KM matrix elements within 
the errors given in Ref. [Q] such that terms contribute constructively. Specifically, for later 
discussions, we use V ud = 0.9741, V us = 0.221, V cd = -0.220, V cs = 0.9740 and V ts = -0.040. 

The charged Higgs boson contribution to the neutron EDM is strongly restricted. The 
dominant term comes from the down quark EDM. Using the valence quark model, we have 

dn ~3 d 3 emdI VaVcFsiml/m^) + T] g V td V ts F 3 (m 2 t /m 2 Hi y 

F ^ x ) = - m X " 5x )( 1 - x ) + ( 4 - 6x ) ln 4 ( 9 ) 

6(1 — x) 6 

Here we have neglected the small QCD correction to the electric dipole operator from the 
gluonic dipole operator induced by operator mixing. The leading QCD correction factor for 
the electric dipole operator is given by |2ljj r/ 7 = [a s (mw) / 'a s (fx)} 16 ^ 33 ~ 2ni '\ 

Using / = (0.35 ~ 3.1) x 10~ 10 GeV~ 2 determined from e'/e, and alowing the light- 
est charged Higgs mass to range from its lower bound around 70 GeV to several hun- 
dred GeV, we estimate the charged Higgs contribution to the neutron EDM as (0.25 ~ 
3.5) x 10 _24 (mrf/300MeV)ecm. Note that d n is proportional to the light quark mass. This 
introduces a further uncertainty because it is not clear whether the current or the constituent 
mass should be used. There are also other uncertainties due to the off shell nature of the 
quarks p5| . 

The neutral Higgs boson exchange gives a contribution which is not well determined. 
Even the sign of the contribution is unknown. Since it is not related to the other parameters 
we have introduced, and since a wide range of the parameter / is still allowed, we will not 
consider the possible contributions from exchange of neutral Higgs bosons in our estimates. 
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VII. B -> 57 IN THE WEINBERG MODEL 



The CP conserving process b — * 57 can place constraints on the CP violating parameters 
of the model P,[lO|, because the CP violating amplitudes contribute to the total rate. In 



the Weinberg model, although dsj interaction is constrained to be small, the corresponding 
657 interaction is enhanced by a factor of (~ m^mb/mlm s )(VtbVts/V cs V c d) ~ 10 5 . Due to 
this enhancement factor, the predicted branching ratio of b — > 37 may be in conflict with 
experimental data. Using the leading log result and normalizing the branching ratio due to 



charged Higgs contribution to the SM one, we have [36 



Br{b -> 57) = 7.1 x KT 4 [(0.313 + 0.273n) 2 + (0.273r 2 ) 2 ], 
WF^ml/m^/Z - Re{a^t)F 2 {ml/w? Hl ) 



ri = l + 



\m{ ai Pt)F 2 {m 2 t /m 2 Hi ) 
F^mljvn^) 



F 1 (x) = — — - — -[(7-5x -8x 2 )(l -x) +x(12- 18s) Ins]. (10) 

\-Zi I -L Ob J 

In the above we have neglected the small contribution from the gluonic bsg interaction. CP 
conserving amplitudes generate the first term in the brackets, and 7*1 contains all of the 
contributions dependent on m t . There are both SM and Weinberg model contributions, and 
there is a region in parameter space where that the CP conserving contributions of the SM 
and of charged Higgs exchange mutually cancel. The CP violating amplitudes generate r 2 , 
which contributes significantly. The branching ratio increases with Higgs mass for fixed /. 
The experimental branching ratio |19|], B(b — > 57) = (3.15 ± 0.54) x 10~ 4 , has recently 



been confirmed by Belle []1T |. For the 95% c.l. upper bound 4.5 x 10~ 4 for b — ► S7 ||19|| . 
we find that there are solutions with |/| < 0.17 x 10~ 10 GeV~ 2 for the charged Higgs mass 
greater than 70 GeV. For larger values of m#, tighter constraints are placed on |/|. Cutting 
the photon energy E 1 to be larger than 2.1 GeV to ensure that the contribution is indeed 
due to the penguin diagram contribution considered here, the central value of the branching 
ratio is reduced slightly to 2.97 x 10~ 4 . The allowed range, at fixed /, for the Higgs mass 
is restricted at the upper end by this reduction. For example, |/| = 0.17 x 10~ 10 GeV -2 is 
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consistent with the reduced branching ratio for 70GeV < van — HOGeV. There is a region 
in which the constraints on / from b — > 37, e and the neutron EDM are consistent. But it 
is not possible to simultaneously satisfy the constraints from b —>■ 57 and from e'/ e. 

VIII. DISCUSSION 

No constraint is yet placed on the model by the results for sin 2/3 from ALEPH, OPAL, 
CDF, BaBar or Belle. 

If the current quark mass m d ~ 10 MeV is used, the resulting value of d n satisfies 
the experimental limit as long as / < 2.56 x 10~ 10 GeV~ 2 . However, if constituent mass 
rrid ~ 300 MeV is used, the model may be in trouble. We know of no convincing argument 
for preferring one mass over the other, and therefore conservatively use the current quark 
mass to estimate limits. 

The values of / from e'/ e, e and d n then have a region of consistency, as do the values 
constrained by B(b — > s-y), e and d n . 

However there is a definite conflict between the limits on / from e'/e and B(b — > 57). 
The latter requires |/| < 0.17 x 10~ 10 GeV~ 2 , and the former requires 0.35 x 10~ 10 GeV~ 2 < 
/ < 3.1 x 10~ 10 GeV~ 2 . As we have been careful to make very conservative estimates of the 
allowed range of / (in the hope of finding that there was still a small region of parameter 
space in which the model is consistent with the data), the gap between these allowed regions 
for / is unbridgeable. Thus we conclude that the Weinberg model is ruled out by the recent 
data for e'/e and for B{b — > 57). 

One of the attractive features of the version of Weinberg model we discuss here is that 
CP violation is generated spontaneously, rather than being put in by hand. If we abandon 
this attractive feature, and explicit CP violation is introduced into the Higgs interaction, as 
in Eq ([[]), as well as into the W interaction by a phase in the KM matrix, this new model 
(which we will call the modified Weinberg model) is not ruled out. The constraint from 
b — > 57 requires that the contribution to e'/e from the Higgs interaction is small, and the 
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main contribution to this CP violating parameter is just the same as that of the SM. As 
has been pointed out in Ref. |]J there are large uncertainties in the SM calculations due 
to our poor understanding of the hadronic matrix elements. One can find allowed regions 
in parameter space in which the experimental value for e'/e is produced in the SM, and 
thus in the modified Weinberg model as well. On the other hand, since the charged Higgs 
exchange can also contribute appreciably to e in this model and can partially cancel the KM 
contribution, as a result the KM phase can have a larger allowed range than that in SM and 
can easily accomodate the larger value of e'/e. Also, in this model, the value of d n can still 
be as large as 3 x 10 _25 (md/300MeV)e cm, very different from the tiny value of the SM [32 
And sin 2/3 can take the large values characteristic of the SM. 
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